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CN ' Abstract. Discrete time evolution of one-dimensional maps is embedded in 

continuous time by truncating the Taylor series expansion of the time evolution 
►> , operator to a finite order N. Truncations with N > 4 leads to unconditional instability. 

(*C) ' Generalization of the truncated models with N = 3 and 4 shows dynamical behaviour 

characteristic of systems with a riddled parameter space 
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1. Introduction 

Time development of systems evolving continuously in time is modelled with differential 
equations, whereas difference equations (maps) are used to describe systems that evolve 
in discrete time. Here we try to study connection, if any, between solutions of discrete 
maps and a set of ordinary differential equations that are derived from the map using 
an ansatz described below. Basically, we truncate the Taylor series expansion of the 
time evolution operator corresponding to the discrete map at a finite order N to get 
an ordinary differential equation (ODE) of order N . Contrary to the expectation that 
the solution of such an ODE should converge to that of the discrete map (when t = an 
integer) in the limit N — > oo we find surprising results. In particular, all these ODEs with 
N > 5 are seen to be unconditionally unstable for almost all initial conditions regardless 
of the details of the underlying map. The nature of the solutions for truncation at N = 3 
and 4 do depend on the details of the map. The 3 rd and 4 th order ODEs obtained for the 
logistic map show period doubling bifurcations leading to chaos, reverse bifurcations, 
instability in some regions etc. Unlike the logistic map, these equations show sudden 
switching from regular to chaos or to unstable behaviour by a slight variation in the 
parameter, which is an indication of structural instability that was discussed in the 
context of systems with a riddled parameter space (Cazelles 2001, Kapitaniak et. al. 
2000, Kim and Lim 2001, Krawiecki and Matyjaskiewicz 2001, Lai and Winslow 1994, 
Lai 2000, Lai and Andrade 2001, Madvinsky et. al. 2001, Maistrenko et. al. 1999, 
Woltering and Markus 2000, Yang 2000, Yang 2001). 

2. Ordinary differential equations corresponding to discrete maps 

Consider a general 1-D map 

x n+1 = f(x n ). (1) 

This discrete time evolution can be embedded in continuous time t by considering the 
equation 

fx(t) = f(x(t)) (2) 

where T is the unit time evolution operator 

t = exp(l) (3) 

We consider the stability of the sequence of equations 

obtained by truncating the formal power series expansion of the time evolution operator 
T to order N. 
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3. Stability analysis 

As mentioned in the introduction, one of the objectives of this paper is to examine the 
nature of solutions of the ordinary differential equations (@J) for various values of N. For 
any map f(x) that is continuous and different iable, we prove that the solutions of 0] are 
linearly unstable for almost all initial conditions if N > 5. 



Consider the stability of solutions of Eq. 
(not necessarily a fixed point) x±. Writing x(t) 

E(T)(^W)+«(z*)^(t> 



j) about an arbitrary reference point 

x* + 5x(t), one gets 
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P{x* 



(5) 



to first order in 5x(t). In Eq. ©, a(x* = 1 - -£f(x)\ x=Xir and f3(x*) = /(x*) - x*. 
Clearly, f3(x+) = if x* is a fixed point of the map. It is easy to show (see appendix A 
for details) that the stability property of the linear inhomogeneous equation (jHJ) is the 
same as that of the homogeneous part 



JV 
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Sx(t) + a(x*)6x(t) = 0. 



J\) [dP 

Since the above equation is linear we look for a solution of the form 
Sx(t) = cexp (fit) 



to get 



3=0 \ J '' 
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We now show that there will exist at least one /i with positive real part irrespective of 
the value of «(£*), if iV > 5. This would imply that the solution of the equation (@J) is 
unstable for any finite N > 5 for almost any initial condition, regardless of the specific 
functional form of f(x). 

Now, Eq. (jHJ) is of the form 

0, (9) 



aofj, + aifi + + oat-i/i + qn 



with aj 
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(N-j)\ 



for < j < (N - 1) and a N 



a(x. 



Nature of zeros of Eq. Q can be examined by using the well known Routh-Hurwitz 
theorem (Korn and Korn 1961). Define 
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Then Routh-Hurwitz theorem states that the number of roots of Eq. (JUJ) with positive 
real parts is equal to the number of sign changes in the sequence {Uj}. If at least one 
sign change occurs in the sequence {Uj}, then Eq. (JOJ) has at least one root (/i) with 
positive real part and hence the solution of Eq. (@J is unstable in the neighbourhood of 
the reference point x*. If this happen for almost any x+, then the solution is unstable 
for any initial condition, globally. 

We now calculate the sequence {Uj} when N > 5. Substituting for a,j in Eq. (JTUJ), 
and on simplifying, we get 

11 2 

Uo= M> ^ 1= (F=1)T' U2 = M(N-2)V (11) 

which are positive for all positive N and 

3 N\(N - 1)\(N - 3)\ y } 

which is negative for all iV > 5. Since there is at least one sign change in the sequence 
{Uj}, it follows that the solutions of Eq. UJ) are linearly unstable for N > 5 irrespective 
of the details of the mapping function f(x), so long as it is continuous and differentiable. 
Since the reference point x* does not explicitly appear in arriving at this conclusion, it 
is evident that the solutions of Eq. (jlj) are unstable for almost any initial condition for 
N > 5. 

We now turn to the case N = 5. Explicit calculation shows that 

TT l TT l TT 2 TT a M Z X ATI («K) - 5/3) 2 + 20/9 

Uo = 5! ' Ul = 4! ' U * = 5!3!' ^ 3 = 5!4! &nd Ua = 5!^ (13) 

From Eq. (|13j) . it is clear that there is a sign change at U3 and no further sign change 
at Ui, if a(x*) < 1. On the otherhand, if a(x ir ) > 1, U3 > and C/4 < 0. Thus there is 
at least one sign change in the sequence {Uj}, regardless of the details of the mapping 
function f(x) for N = 4, as well. No such general conclusions can be drawn for the 
cases N = 1, 2, 3 and 4. 

We now turn to the specific case of truncation of chaotic maps. Since a minimum 
of three dimensions is required for occurrence of chaos in ordinary differential equations, 
truncations to order N — 1 and 2 are uninteresting from this perspective. We also now 
know that the truncations to order N > 5 lead to unconditional instability. Therefore, 
in what follows, we concentrate on truncations with N = 3 and 4. It is found (Rajan 
Nambiar 2003) that the solutions of these equations show many features which are not 
shared by the solutions of the original discrete map. For concreteness, we present the 
results for the logistic map f(x) = px(l — x), x e [0,1]. 

4. Truncation at N=3 and 4 

For this case, Eq. (J1J) becomes 



Stability of differential equations associated with maps 5 

Linear stability analysis of the above equation, with f(x) = px(l — x) shows that both 
fixed points (0 and 1 — 1/p) are stable for p < 4. Thus the truncation at N = 3 leads 
to regular behaviour for all p e [0,4], the range of p for which the logistic map is map 
of the unit interval to itself. This has to be contrasted with the occurrence of chaos in 
the logistic map when p e [3.66..., 4]. However, the system of equations (J14j) do show 
chaotic behaviour when p > 4. 

Using the scaled variables X— (2p/9)x and r =t/3, Eq. (fT4^) can be rewritten as 
afX d 2 X dX , „ 

^ + ^ + "* + - xx + x =° < 15) 

where v = 2/3 and A = 2{p — l)/9. In fact, the generalization of Eq. (fTHjl obtained by 
allowing v to take arbitrary values is equivalent to the equations studied by Coulett et. 
al. (Coulett et. al. 1979) and Arneodo et. al. (Arneodo et. al. 1985) which exhibited 
striking features in its solution. The model shows regular, chaotic and or unstable 
behaviour for certain choices of the parameters v and A. For example, for a fixed z/, as A 
is increased, we observe sequences of either finite or infinite number of period doubling 
bifurcations and reverse bifurcations with unstable regions in between. The bifurcation 
diagram is substantially different for a neighbouring v. Thus the significant aspect of 
this equation is that the nature of solution changes drastically for very small changes in 
the values of the parameters. It is as if the parameter space is riddled. The solution of 
ODE obtained by truncating at N = 4 also show similar features. 

5. Conclusions 

The present work shows that the nature of solutions of the ordinary differential 
equations, obtained by truncating the power series of the time evolution operator 
corresponding to discrete maps, are very different from those of the discrete maps. 
In particular, the solutions with N > 5 are unconditionally unstable regardless of the 
details of the map. For the specific case of the logistic map, truncations at N = 3 and 
4 lead to solutions characteristic of systems with a riddled parameter space. 
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Appendix A. Solution of the inhomogeneous equation (HjJ) 



Let us define an N- dimensional column vector £ 
£ = (fi, -,6v) T , with & = Sx(t), £ i+ i = — £ i; j 

Equation © can then be rewritten as 
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The formal solution of Eq. (JA.2J) is given by 

f(t) = exp [Mt]to + / dt'exp [M(t - t')]v{t') 
Jo 

Using the fact that v is independent of time and redefining the variable of integration, 
we get 

(A.5) 
(A.6) 



£(t) = exp [M£]f + ( / dt'exp [Mf] J v 

Let S be the similarity transformation matrix that diagonalises M so that 

S _1 MS = D = Diag(/j,i, , //at) 

It then follows that 

(pMl* _ 1 pV-Nt _ i \ 
, , S-V (A.7) 
\X\ [1 N J 

The first term on the right hand side of Eq. (|A.7|) corresponds to the solution of the 
homogeneous part (see Eq. ©)• It is clear that if at least one of the {//,■} has a positive 
real part so that the solution of the homogeneous equation is unstable, so is the solution 
of the inhomogeneous equation. 



